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Recently proposed scenarios for the cuprates make extensive use of a "flat" quasiparticle (q.p.) 
dispersion and short-range hole-hole interactions in real-space, both caused by antiferromagnetic 
(AF) correlations. The density of states (DOS) at half-filling has a robust peak which boosts the 
superconducting critical temperature T c to large values as holes are introduced into the (rigid) 
q.p. band. Here, the stability of such scenarios is studied after a finite but small hole density 
is introduced. The overall conclusion is that the main features of real-space AF-based approaches 
remain qualitatively similar, namely a large T c is found and superconductivity (SC) appears in the 
d x 2_ y 2 channel. As the hole density grows the chemical potential /j, crosses a broad peak in the 
DOS. We also observe that extended s-wave SC competes with d-wave in the overdoped regime. 
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The presence of "flat" regions in the normal-state q.p. 
dispersion is a remarkable feature of the phenomenol- 
ogy of hole-doped high temperature superconductors. 0] 
These flat bands are located around momenta p = (tt, 0) 
and (0, 7r), and at optimal doping they are ~ lOmeV be- 
low the Fermi energy, according to angle-resolved pho- 
toemission (ARPES) studies, jl] Antiferromagnetic cor- 
relations may play an important role in the generation of 
these features, as has been suggested by studies of holes 
in the 2D t- J and Hubbard models at and away from half- 
filling. In addition, the flat regions induce a large 
peak in the DOS which can be used to boost T c . This 
leads to a natural explanation for the existence of an "op- 
timal doping" which in this framework occurs when the 
flat regions are reached by /i. We will refer to these ideas 
as the "Antiferromagnetic van Hove" (AFVH) scenario. 
H Superconductivity in the d x 2„ y 2 channel is natural in 
the AFVH scenario due to the strong AF correlations. 
[0,0 The interaction of holes is better visualized in real 
space i.e. with pairing occuring when dressed holes share 
a spin polaronic cloud, as in the spin-bag mechanism. ||] 
This real-space picture (see also Ref. 0j) holds even for 
a small AF correlation length, S,af, and in this approach 
there is no need to tune parameters to work very close to 
an AF instability. Previous scenarios have also used van 
Hove (yH) singularities in the band structure to increase 
T c , jl0[ but d-wave SC is not natural in this context un- 
less AF correlations are included. 

To obtain quantitative information from these intuitive 
ideas, holes moving with a dispersion calculated using 
one hole in an AF background, caf(jp), and interacting 
through a nearest-neighbors (NN) attractive potential, 
that mimics the sharing of spin polarons, have been pre- 
viously analyzed. || Within a rigid band filling of eaf(p) 
and using a BCS formalism, d x 2_ y 2 superconductivity 
dominates with T c ~ 100K caused by the large DOS im- 
plicit in the hole dispersion. However, the accuracy of 
the rigid band approximation has not been studied be- 
fore. In particular, the following questions arise: (i) does 



the q.p. peak in the DOS found at (n) = 1 survive a 
finite hole density?; (ii) how does the q.p. dispersion 
change with hole doping, and to what extend does this 
affect previous calculations in this framework?; (iii) are 
the "shadow" regions generated by AF correlations [p|Jll|l 
important for real-space pairing approaches? In this pa- 
per all these issues are discussed. The overall conclusion 
is that as long as the hole density is such that the £af is 
at least of a couple of lattice spacings, the predictions of 
the original AFVH scenario || and other similar theories 
remain qualitatively the same. 

To study the presence of a large DOS peak away from 
half-filling, N(u>)(= J2 P ^(p> w )) nas been calculated nu- 
merically in systems with strong AF correlations. While 
accurate extrapolations to the bulk limit are difficult, the 
simple qualitative picture emerging from these studies 
seems robust. In Fig. la, N(uj) for the 2D t-J model 
obtained with exact diagonalization (ED) techniques is 
shown at several densities. |?J At half-filling, a large DOS 
peak caused by flat regions in the q.p. dispersion appears 
at the top of the valence band, as discussed before. ||] 
Weight exists at energies far from i-e. the large peak 
carries only part of the total weight. As shown in Ref. 
the maximum in the DOS is not reached at the top 
of the valence band but at slightly smaller energies. As 
(n) decreases, the peak is now much broader but it re- 
mains well-defined. At (n) = 0.87, [i is located close to 
the energy where N(u>) is maximized. If a source of hole 
attraction exists in this system, a SC gap would open at 
fx, and the resulting T c could be large. However, even 
with T c ~ lOOi-T the SC gap would not alter much the 
already robust peak in the DOS shown in Fig. la. For 
(n) ~ 0.75, fi moves to the left of the peak, where T c 
should become smaller. Then, as (n) is reduced away 
from half- filling in the 2D t-J model Fig. la suggests that 
/j, travels across a broad DOS peak providing a natural 
definition for the "undcrdoped, optimal and overdoped" 
densities. [|l2| Note that since the peak width increases 
substantially with hole density, strictly speaking the rigid 
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band filling is invalid. However, such an approximation 
may have captured the qualitative physics of the prob- 
lem, as discussed in this paper, since the DOS peak is 
not washed out by a finite hole density. 
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FIG. 1. (a) N(uj) for the 2D t-J model averaging results for 
clusters with 16 and 18 sites, at J/t = 0.4 and for the densi- 
ties indicated. The 5-functions were given a width rj — 0.25i. 
Similar results were found at other values of J/t; (b) N(u) 
for the one band Hubbard model obtained on a 4 x 4 clus- 
ter using QMC and ME. The temperature is T = t/4 and 
U/t = 12. Densities are indicated. Previous studies showed 
that QMC/ME results are only weakly affected by size ef- 
fects. [3] 

Quantum Monte Carlo (QMC) data for the 2D one- 
band Hubbard model in strong coupling indicate a qual- 
itative behavior similar to that of the t-J model, i.e. u 
also crosses a peak as (n) varies. This effect was observed 
in previous simulations, |l3f| and also in our own studies 
(Fig. lb). |Q However, note that QMC simulations must 
be supplemented by a Maximum Entropy (ME) analysis 
to obtain N(u>). ME predictions are more accurate close 
to bj than at finite frequency, and thus we here do 
not interpret Fig. lb as suggesting the "generation" by 
doping of a Kondo-like peak at the top of the valence 
band which does not exist at (n) = 1. We believe that 
the absence of a sharp peak at half-filling is caused by 
the systematic errors of the ME procedure and the fi- 
nite temperature of the simulation. Actually, a robust 
peak at (n) = 1 appears in other computational studies, 
P,[l3| and in the t-J model (Fig. la). Experimentally in 
the cuprates it has been already established that the 
states observed in PES upon doping are already present 
in the insulator and are not Kondo resonances. In the 
present calculation, Fig. la suggests that the broad peak 
observed at (n) < 1 in the t-J model is smoothly con- 
nected to the sharper peak found at half-filling. 

Similar conclusions can be reached in geometries other 
than the 2D square lattice. In particular, we studied the 
DOS of the t-J model on a ladder, finding results very 



similar to those of Fig. la, in agreement with previous cal- 
culations. (l6) Note that on ladders the presence of pair- 
ing at finite hole density has been predicted. Q Thus, 
a gap may open at ii once couplings and lattice sizes are 
reached where pairing effects are important. We have ob- 
served such effects in our studies, but we also found that 
at (n) = 1 the DOS has a sharp peak, as in Fig. la. Then, 
the ladder DOS peak is not exclusively caused by a BCS- 
induced redistribution of spectral weight upon doping. 

In all cases discussed here, calculations of the spin-spin 
correlations show that ^af is approximately a couple of 
lattice spacings when \i is located near the DOS peak, 
becoming negligible as the overdoped regime is reached. 
Thus, a nonzero £af and \i near a large DOS peak are 
correlated features. It is in this respect that scenarios 
where AF correlations produce a large DOS that en- 
hances T c differ from vH theories where divergences in 
the DOS are caused by band effects already present be- 
fore interactions are switch on. 

While the existence of a robust peak in the DOS is in 
good agreement with ARPES data, 0] it is in apparent 
disagreement with specific heat studies for YBCO which 
have been interpreted as corresponding to a flat DOS. 
p8| The lack of p-resolution in the specific heat measure- 
ments may solve this puzzle. Actually, angle-integrated 
PES results for Bi2212 do not show the sharp flat features 
found in ARPES for the same material. |l5|JT9| ] Similar 
effects may affect the specific heat data. 

Now let us discuss the p-dependence of the q.p. band 
obtained from A(p, u>). Results are already available in 
the literature, pjffiifl In Fig.2a-b, the q.p. dispersion is 
shown at (n) — 1, and at (n) ~ 0.87 for the t-J model. 

Note that upon doping the small q.p. bandwidth and 
the flat regions remain well-defined, inducing a large DOS 
peak (Figs. la). However, the region around p = (7T,7r) 
has changed substantially, i.e. the AF shadow region 
clearly observed at (n) — 1 has reduced its intensity and 
considerable weight has been transferred to the inverse 
PES region (Fig. 2b). Actually, the q.p. dispersion at 
(n) ~ 0.87 can be fit by a tight-binding nearest-neighbors 
(NN) dispersion with a small effective hopping likely as- 
sociated with J. This effect was also noticed in studies 
of the Hubbard model. |l|J|f§. 

The changes in the q.p. dispersion with hole doping 
can be interpreted in two ways: First, note that A(p, ui) 
is influenced by matrix elements of bare fermionic opera- 
tors connecting states with N and N ± 1 particles. This 
is important when the q.p. weight is small i.e. when 
the state c pa \gs)N does not have a large overlap with 
the ground state of the N — 1 particles subspace \gs) n-i 
(\gs)N being the ground state with N particles). It has 
been suggested |^,^3| that if the hole excitation is cre- 
ated by a new operator 7 po - that incorporates the dress- 
ing of the hole by spin fluctuations, then r y pa -\gs)N may 
now have a large overlap with \gs)N—i- In other words, 
if the dressed hole state resembles an extended spin po- 
laron, then the physics deduced from PES studies, which 
rely on the sudden removal of a bare electron from the 
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system, may be misleading. To the extent that spin po- 
larons remain well-defined at finite density, the use of 7 pcr 
will induce spectral weight rearrangements between the 
PES and IPES regions, and Fig. 2b could resemble Fig. 2a 
after such spectral weight redistribution takes place, as 
some results already indicate. [^3| If this idea were cor- 
rect, then the study of SC and transport, both regulated 
by dressed quasiparticles, could indeed be handled by 
filling a rigid band given by €af(p)- 
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is justified. Note also that the local character of the real- 
space potential does not imply small Cooper pairs. Their 
size can be regulated using both its range and intensity. 
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FIG. 2. (a) q.p. energy vs momentum obtained at 
half-filling. [2] The result shown, 6af(p), is a good fit of Monte 
Carlo data on a 12 x 12 cluster at J/t = 0.4; (b) q.p. disper- 
sion vs momentum at (n) = 0.87 and J/t = 0.4 using exact 
diagonalization of 16 and 18 sites clusters (from Ref. [19]). 
The open (full) circles are IPES (PES) results. Their size is 
proportional to the peak intensity. The solid line is the fit 
£]vjv(p); (d) V(x) along the x-axis after Fourier transforming 
the smeared potential V(p) = 5(p — Q) (see text), (af is 
given in lattice units. 

However, an alternative is that the difference between 
Figs. 2a and 2b corresponds to an intrinsic change in 
the q.p. dispersion as ^af decreases. This is the less 
favorable case for the AFVH approach, and thus we 
analyze this possibility in detail here. For this pur- 
pose, we have applied the standard BCS formalism to 
a model with a low density of q.p.'s having a disper- 
sion 67VAr(p)/eP = — 0.2(cospa; + cosp y ) , which roughly 
reproduces the dominant features of Fig. 2b. As before, 
we include NN attraction between q.p.'s. While naively 
it may seem dubious to use the same interaction both 
at and away from half-filling, the hole-hole potential 
should not be much affected at distances shorter than 
£af< This can be illustrated by a real space analysis 
(Fig. 2c) of a smeared (^-function potential of AF origin 
V(q) = + £ AF (q — Q) 2 ], where the lattice spacing 

is set to 1, and Q = (ir, tt). Fig. 2c shows that the NN 
potential (x = 1) does not change noticeably as £af is 
reduced, while V(x > 1) is rapidly suppressed. Then, us- 
ing the same NN form of the potential for many densities 
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FIG. 3. (a) T c vs u obtained with the BCS gap equation, 
with u = as the chemical potential corresponding to the 
saddle-point. The dashed line corresponds to eaf(p), and 
the solid line to ejvjv(p). Interpolations between these two 
extreme cases can be easily constructed using two band mod- 
els with weights regulated by Z-factors. Note the presence of 
both d x 2_ y 2-wave and extended s-wave SC; (b) T c for d-wave 
SC vs the percentage x of filling of the q.p. band (not of the 
full hole spectrum). The dashed line is the same as in (a). 
The solid line corresponds to the BCS gap equation result 
making zero the weight Z p of states in the q.p. dispersion 
that are at energies from the saddle point larger than 2.5% 
of the total bandwidth (i.e. basically including states only in 
a window of energy ~ 125K around the flat regions). The 
dot-dashed line is the same but using a window of ~ 250K 
around the flat regions. 

Then, let us discuss the results obtained using the BCS 
formalism with an attractive NN potential V = — 0.6J, 
II with J = 0.125eP, for both dispersions caf(p) and 
ejvw(p)- Solving numerically the gap equation, T c is 
shown in Fig. 3a. As reported before, || SC at T c ~ 
80 — 100-ftT in the d x 2_ y 2 channel appears naturally if 
the AF-dispersion is used. If, instead, ejv n (p) is used, 
the flat bands present in this narrow dispersion produce 
T c ~ 30 K which is still large. Even more remarkable 
is the fact that the d x 2_ y 2 character of the SC state is 
maintained. This result can be understood noticing that 
a combination of £nn(p) with an attractive NN poten- 
tial effectively locates us in the family of "t-U-V" models 
with U repulsive and V attractive, where it is known 
that for a "half-filled" electronic band the dominant SC 
state is d x 2_ y 2-wave. [p4| In other words, when /i is at 
the flat region in Fig. 2b it approximately corresponds to 
a "half-filled" q.p. band, leading to a d x 2_ y 2-wave SC 
state (Fig. 3a). |25| Then, one of our main results is that 
even if the q.p. dispersion changes substantially with 
doping near the Q point, such an effect does not alter 
the main qualitative features found in previous studies. 
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The present analysis also shows that the AF "shadow" 
regions of £af(p) are not crucial for the success of the 
real-space approach. Actually, using 6_/v./v (p) j which docs 
not contain weight in PES near T c is still robust 

and the d-wave state remains stable. To make this point 
more clear, we analyzed T c using €af(p) but modulating 
the contribution of each momentum with a p-dependent 
weight Z p in the one particle Green's function. We con- 
sidered the special case where Z p is zero away from a win- 
dow of total width W centered at the saddle point, which 
is located in the flat bands region. Inside the window W, 
Z p = 1. Such a calculation also addresses indirectly pos- 
sible concerns associated with widths ~ (e(p) — e^) 2 that 
q.p. peaks would acquire away from half-filling in stan- 
dard Fermi liquids. Results are shown in Fig. 3b, for d- 
wave SC. Note that even in the case where W is as small 
as just 5% of the total bandwidth (itself already small of 
order 2J), T c remains robust and close to 70K. Then, it 
is clear that the dominant contribution to T c comes from 
the flat regions and the shape of the q.p. dispersion away 
from them is not relevant for the sucess of the real-space 
approach. 

We end this paper with novel predictions obtained 
from 6af{p) and 6nn(p), when a reaches the bottom 
of these bands (which should not be confused with the 
bottom of the full hole spectrum since a large amount of 
weight lies at energies lower than those of the q.p. band). 
In this regime, the BCS analysis applied to any of the 
two dispersions shows that extended s-wave SC domi- 
nates over d x 2_ y 2-wave SC in the "overdoped" regime 
(Fig. 3a). This occurs when the q.p. band is nearly empty 
which, according to Fig. la, corresponds to an overall den- 
sity of (n) ~ 0.7. This change in the symmetry of the 
SC state can be understood recalling again that €nn(p)- 
supplemented by a NN attraction formally corresponds 
to an effective "t-U-V" model. It is well-known that in 
this model the SC state symmetry changes from d- to 
s-wave as the density is reduced away from half-filling 
to a nearly empty system. ^J,^6| Actually the bound 
state of two particles with a NN tight binding disper- 
sion and NN attraction is s — wave. To the extend that 
the AF- or NN-dispersions survive up to ~ 25% hole 
doping, as suggested by numerical data, (l3|,^0| scenar- 
ios based on the real-space interaction of q.p.'s predict 
a competition between extended s-wave and d x 2_ y 2-wave 
SC in the overdoped regime. It is remarkable that calcu- 
lations based on the analysis of the low electronic density 
(n) <C 1 limit of the t-J model lead to analogous conclu- 
sions. 26 Our approach is based on a very different 
formalism, but it arrives to similar results, and thus we 
believe that a crossover from d- to s-wave dominated SC 
in overdoped cuprates could occur. Recent ARPES data 
for overdoped Bi2212 have indeed been interpreted as 
corresponding to a mixing of s- and d-wave states. p8| 

Summarizing, the predictions of real-space theories 
based on dispersions and interactions calculated at half- 
filling arc qualitatively stable upon the introduction of 



a finite hole density. These include a large T c , within 
the BCS formalism, and SC in the d x 2_ y 2 channel. In 
the overdoped regime a possible crossover from d- to ex- 
tended s-wave SC was discussed. 
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